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We present Monte-Carlo simulations of jet evolution in lead-lead collisions at the Large Hadron Collider (LHC) 
at CERN focusing on the dijet asymmetry measured by the ATLAS and CMS collaborations. In the simulation, 
T~H hard partons are interacting with the hydrodynamical background medium, undergoing radiative and coUisional 

C 2 ^ processes. The measured dijet asymmetry is well described by the simulation and is hence consistent with 

^N) partonic energy loss in a hot, strongly-interacting medium. 

Oh 
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^ 1. Introduction 

1^ 

With a center of mass energy of -y/s]vjv = 2.76 TeV at the LHC, which is significantly larger than that achieved 
^ at the RHIC, far more energetic jets are kinematically accessible at the LHC. The ATLAS collaboration was 

(~| able to measure over 1000 dijets where the leading jet has transverse energy Et > 100 GeV and the opposing 
Qh jet has energy Et > 25 CeV [1]. The CMS collaboration performed a similar analysis on a large sample of jets 
{Eti > 120 GeV, Et2 > 50 CeV) [5]. These results are a significant improvement over the results from the 
RHIC, where the total energies of the jets were far lower and therefore harder to separate from fluctuations in 
the underlying bulk. Also, the models for partonic evolution rely on the probe parton having high energy, and 
when this separation of energies exists one can expect the hadronization of these partons to be described well 
with vacuum fragmentation functions. 
^ Several theoretical studies have addressed the measured dijet asymmetry at the LHC. In 15] the evolution of 

CN| the jet shower was studied, analyzing the jet propagation through the quark-gluon plasma and its interaction 
^\ with the medium. Cood agreement with the experimental data was found. The authors of [1] conclude that the 
removal of soft components from within the jet cone via elastic collisions will induce a dijet asymmetry. Also, 
the PYQUEN model [5j was used the to quantify the "jet-trimming" . 
On In this work, we present a somewhat more detailed analysis of what was presented in |6], where we applied 

MARTINI to lead- lead collisions at the LHC [3 [5] . In Section [2j the physics behind martini is reviewed, as 
well as the description of the bulk of heavy-ion collisions with 3+1-dimensional hydrodynamics. In Section [3j 
. . runs of MARTINI with cuts given by the ATLAS and CMS detectors and their analyses are compared with the 
experimental results for both dN/dAj (the yield of dijets differential in v4j, where Aj measures the energy 
anisotropy of dijets) and dN/dcj). 
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2. Transport of high-energy partons and martini 

MARTINI solves the rate equations 
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where the various differential rates dT\^^{p^k) / dkdt determine the spUtting of partons I and m, one with mo- 
mentum fc, from a parton i with momentum p 9J. 

In its current implementation, martini uses rates which take into account both radiative and cohisional QCD 
processes, calculated at finite temperature. Collisional processes involve soft momentum transfers sensitive to 
the gluon's screening mass and therefore, hard thermal loop results at leading order are used to describe these 
processes. For the elastic processes martini does not depend on the "diffusive approximation" : there is no need 
to assume that the rates are only significant when uj is small p[0 . It is important to point out at this point 
that for the elastic collisions there is no exact energy and momentum conservation within the hard partons 
because energy is lost to the strongly coupled hydrodynamic medium. Because at this point this lost energy 
and momentum is not added to the energy in the hydrodynamic calculation (the hydrodynamic background is 
not modified by the hard partons), it is lost from the analysis. The underlying assumption here is that this 
energy and momentum is thermalized rapidly and becomes part of the background medium. Radiative energy 
loss is modeled using the Arnold-Moore- Yaffe approach, where the interference of bremsstrahlung gluons from 
multiple scatterings is taken into account with an LPM-like integral equation for the energy loss rate [TTUl3j . 
Also for the radiative processes a certain cut in momentum has to be included, because the AMY formalism 
is not applicable for parton energies of the order of the temperature. We assume that radiated partons with 
momenta below 2 GeV are better described as part of the strongly coupled medium and do not integrate them 
into the hard parton shower. Furthermore, partons with energies below 4T do not undergo any perturbative 
interactions anymore. A more detailed analysis of the dependence of final results on these separation scales 
between the perturbative regime and the strongly coupled bulk medium is currently on the way and will be 
presented in a forthcoming work. 

The momenta of high-energy partons are sampled using pythia event generation |14| . and their initial 
positions in the transverse plane of heavy-ion collisions are sampled according to ns(x, y, b), the distribution of 
binary collisions for a given impact parameter b of the collision. These partons are then evolved through the 
background of bulk particles. In this work, this evolving background is modeled using MUSIC, a 3+1-dimensional 
hydrodynamic simulation |15j . 

For the results in Section [3j martini is run with Us — 0.27 including both collisional and radiative processes, 
as well as with radiative and collisional energy loss separately. The finite-temperature rates for these processes 
are determined by the temperatures and flow in lead-lead collisions as simulated with MUSIC for an impact 
parameter of 6 = 2.31 fm, reproducing the multiplicities of the 0-10% centrality class. In this study we use a 
simulation with ideal hydrodynamics starting with averaged initial conditions. 



3. Results for lead-lead collisions measured at ATLAS and CMS 

Once high-energy partons have evolved and hadronized, the resulting hadrons must then be reconstructed into 
jets. For the best possible comparison with the results of the LHC, we use the same anti-fct jet reconstruction 
that the ATLAS collaboration uses [TB]. These algorithms depend on the definition of distances between two 
4- moment a: 



The distances are determined between all pairs of final-state particles whose energies are large enough to trigger 
the calorimeters, and starting with the smallest distance, 4-momenta close to each other are clustered and added 
together and final jets are determined. The implementation of this algorithm that we used is fastjet, publicly 
available online [17) . 

Once the clustering of hadrons into jets is complete, the jet with highest is determined, and the highest 
energy jet whose azimuthal angle from the leading jet Acj) > 7r/2 (or 27r/3, as is the case with the CMS analysis) 
is also determined. If the energies of this dijet are high enough to make it into the given detector's analysis, 
they are recorded and binned. 

In Figure[l] we show the results for ATLAS, in the 0-10% centrality range, for the differential yield dN/dAj, 
where Aj — ^ti-Et2 measure of the transverse energy asymmetry of the dijets. The ATLAS results 
used are from the latest analysis using R = 0.4 [TSj; there was little dependence on R found in the latest 
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Figure 1: The difFerential yield dN/dAj for proton-proton collisions and lead-lead collisions. 



results, suggesting partonic energy loss as the dominant mechanism leading to dijet asymmetry. Our results 
are compared with p-l-p events using pythia and fastjet, and the differential yields are normalized to one. In 
Figure [2j we show the difFerential yields dN/dip, where is the azimuthal opening angle for the dijets. 

The presented results are for as — 0.27. Fig. [2] shows no significant difference in the distribution of dijets 
between proton-proton and lead-lead collisions. The experimental results show an increase in the yield at small 4> 
in lead-lead collisions over what was observed in proton-proton collisions. This enhancement, while significant, 
affects a relatively small number of dijets in ATLAS' sample, and can be explained by uncertainties in the 
combinatorics. For example, if one of the jets is absorbed and an uncorrelated background fluctuation (or part 
of a second dijet) is used as associated jet instead, this could lead to such enhancement of dN/dip. Of course 
it will be most noticeable at Acj) ^ 7r/2, where dN/dcfi is smaller. To get a feeling for the relevant ingredients 
needed to reproduce the experimentally observed dijet asymmetry, we present results of calculations with only 
elastic and only radiative processes in Figs. |3]and[4j respectively. 

The reason why there is a difference to the proton-proton result in both of those cases is the following. Partons 
can get kicked out of the cone due to elastic collisions and hadrons can get kicked out due to the fragmentation 
fcr in the first case. There is also some additional asymmetry because the energy of the shower is not completely 
conserved in the calculation - energy and momentum is lost to the bulk medium and we do not keep track of it, 
assuming that it thermalizes in the strongly coupled medium. In the case of only radiative processes, there is 
broadening due to the fragmentation kx, which can kick hadrons coming from soft radiated partons out of the 
cone, but also the previously mentioned cutoff below which radiated gluons are not added to the shower causes 
some asymmetry. A detailed study of the dependence on this cutoff is on the way. 

Together, this leads to a description of the process in the full simulation. It is apparent from the figures 
that both radiative and elastic processes are needed to describe the experimentally observed dijet asymmetry. 
Radiative processes generate soft partons which can be kicked out of the cone more easily than hard partons by 
elastic collisions. Furthermore, the fragmentation fc-p in pythia's Lund model will also transport hadrons from 
softer partons out of the cone more efficiently. 

The soft physics of the shower partons is not described in detail within this pQCD based simulation such that 
some assumptions, manifest in the low momentum cutoffs for certain processes, had to be made. Nevertheless, 
the general physical process can be well reproduced with this perturbative description of the hard degrees of 
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Figure 2: The differential yield dN/dAcfi for proton-proton collisions and lead-lead collisions. 
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Figure 3: The differential yield dN/dAj for proton-proton collisions and lead-lead collisions, including only elastic 
processes. 
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Figure 4: The difFerential yield dN/dAj for proton-proton collisions and lead-lead collisions, including only radiative 
processes. 



freedom coupled to a strongly interacting bulk medium. 

Figure [5] shows the differential yield in Aj determined by CMS' dijet sample, compared with martini's results 
based on C!MS' kinematical cuts 



4. Conclusions 



The study reported here utilizes the pQCD and thermal-QCD based MARTINI numerical simulation with a 
hydrodynamic background determined by MUSIC and full jet reconstruction using fastjet. Using only one free 
parameter - as - we can explain a large part of the jet asymmetries observed in the recent ATLAS and CMS 
experiments at the LHC as the consequences of high energy jets interacting with the evolving QGP medium. 
The alignment of the dijets dN/d(p is also well reproduced up to effects resulting from the fluctuating background 
that we do not yet include. Further details and computation of the jet fragmentation functions will be presented 
in a future work. 
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Figure 5: The differential yield dN/dAj for lead- lead collisions at ^/s — 2.76 GeV compared to CMS data [2]. 
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